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ABSTRACT 

The full text of the discussion paper presented at the Whither Turbulence Work- 
shop (Cornell University, March 22-24, 1989) on past and future trends in turbulence 
modeling is provided. It is argued that Reynolds stress models are likely to remain the 
preferred approach for technological applications for at least the next few decades. In 
general agreement with the Launder position paper, it is further argued that among 
the variety of Reynolds stress models in use, second-order closures constitute by far 
the most promising approach. However, some needed improvements in the specifica- 
tion of the turbulent length scale are emphasized. The central points of the paper 
are illustrated by examples from homogeneous turbulence. 
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Discussion of Turbulence Modeling: Past 
and Future 

Charles G. Speziale 

Institute for Computer Applications in Science and Engineering 

NASA Langley Research Center 
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1. Introductory Remarks 

The thrust of the position paper by Launder is that second-order closure models 
represent the best hope for the reliable prediction of the complex turbulent flows of 
technological interest both now and in the foreseeable future. By building on the pio- 
neering research of Rotta [1] and by introducing some fundamental new ideas, the work 
of Launder, Lumley, and others has without doubt made significant contributions to the 
advancement of second-order closures. In the position paper by Launder, a strong case 
is made for the superior predictive capabilities of second-order closures in comparison to 
two-equation models or eddy viscosity models. Most notably, turbulent flows involving 
rotations and streamline curvature have been shown by Launder and others [2,3] to be 
better described by second-order closure models. The same is true for turbulent flows 
with stratification and relaxation effects. Launder very aptly cites four active areas of 
research for the improvement of second-order closures: (i) models for the rapid pres- 
sure strain correlation, (ii) models for the turbulent diffusion terms, (iii) adjustments 
for near wall turbulence effects, and (iv) modeled transport equations for the turbulent 
dissipation rate or length scale. 

In the sections to follow, by making use of some simple examples from homogeneous 
turbulence, the primary point made by Launder concerning the superior predictive ca- 
pabilities of second-order closures will be amplified. Most notably, it will be shown that 
second-order closures are capable of describing the stabilizing or destabilizing effect of 
rotations on shear flow - a problem which cannot be even remotely analyzed by the 
simpler models. However, some lingering problems concerning the development of ad- 
equate models for the rapid pressure-strain correlation and the turbulent length scale 
will be emphasized (see Speziale [4]). In regard to the latter issue, the strengths and 
weaknesses of the commonly used modeled dissipation rate transport equation will be 
discussed and a definitive argument will be put forth as to why previous attempts at the 
development of an improved dissipation rate transport equation have failed. Alternative 
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approaches based on a tensor length scale will also be discussed along with the author’s 
views concerning the prospects for future research. 

2. The Case for Second- Order Closure Models 

The commonly used eddy viscosity models and two-equation models have three major 
deficiencies [5,6]: 

(i) the inability to properly account for rotational strains, 

(ii) the inaccurate prediction of normal Reynolds stress anisotropies, 

(iii) the inability to account for component Reynolds stress relaxation 
and amplification effects. 

In so far as point (i) is concerned, it should be noted that the K- e model is oblivious 
to the presence of rotational strains (e.g., it fails to distinguish between the physically 
distinct cases of plane strain, plane shear, and rotating plane shear). Other commonly 
used algebraic eddy viscosity models such as the Baldwin- Lomax model are also funda- 
mentally incapable of describing the effect of rotations on sheared or strained turbulent 
flows [2,6]. As alluded to in point (ii), all eddy viscosity models, including the K-e 
model, yield highly inaccurate predictions for the normal Reynolds stress anisotropies in 
simple turbulent shear flows. This makes it impossible to describe a variety of secondary 
flow phenomena (e.g., the K-e model erroneously predicts that there are unidirectional 
mean turbulent flows in non- circular ducts in contradiction to experiments which indi- 
cate the presence of an additional secondary flow [5]; see Figure 1). These problems can 
be partially overcome by the use of two-equation turbulence models with a nonlinear 
algebraic Reynolds stress model (see Launder and Ying [7], Rodi [8], and Speziale [5]), 
but only for turbulent flows that are nearly in equilibrium. Non-equilibrium turbulent 
flows that have a spatially or temporally evolving structure (e.g., the flows with relax- 
ation or amplification effects mentioned in point (iii)) cannot, in general, be described 
properly by two-equation models. For example, in an initially anisotropic turbulence, 
where at some time t = t 0 the mean velocity gradients are set to zero, the K-e model 
erroneously predicts an instantaneous return to isotropy wherein 

Tij = KSij , t > t 0 (1) 

or equivalently, 

kj = 0, t>t 0 (2) 

given that K is the turbulent kinetic energy, = — u^uj is the Reynolds stress tensor, 
and bij = — (ry + jK8ij)/2K is the anisotropy tensor. In considerable contradiction 
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to (1), experiments indicate that there is a very gradual return to isotropy - an effect 
that can be characterized much better by second-order closure models. In Figure 2, 
the temporal evolution of the second invariant of the anisotropy tensor II is shown 
corresponding to a relaxation from the plane strain experiment of Choi and Lumley 
[9]. From this graph, it is clear that the second-order closure model (which is based on 
the Rotta model for the slow pressure-strain correlation) does a reasonably good job in 
reproducing the experimental trends unlike the K-e model which erroneously predicts 
that II = 0 for dimensionless time r > 0. 

Now, the greater predictive capabilities of second-order closure models will be demon- 
strated by a simple, but non-trivial, example which is not often discussed in the turbu- 
lence modeling literature. The problem to be considered is homogeneous turbulent shear 
flow in a rotating frame (see Figure 3). This problem constitutes a non-trivial test of 
turbulence models since it involves arbitrary combinations of shear and rotation which 
can have either a stabilizing or destabilizing effect. 

For any homogeneous turbulent flow, the standard K-e model takes the general form 


[ 10 ] 


r * “ ~Z K8ij + Cfl e ( dx, + dx? 
• dv ; 

k = t '^- c 

6 &Ui £ 2 

*- c «7»l£- c *K 


( 3 ) 

( 4 ) 

( 5 ) 


in all frames of references independent of whether or not they are inertial. In (3)- 
(5), V{ is the mean velocity field, e is the turbulent dissipation rate, and C^, C ei and 
C '2 are constants which assume the values of 0.09, 1.44, and 1.92, respectively. For 
homogeneous turbulent shear flow in a rotating frame (as shown in Figure 3), the mean 
velocity gradient tensor is given by 


dv{ 

dxj 


0 S 0 \ 
0 0 0 
0 0 0 / 


( 6 ) 


and 12 { = (0, 0, Q) is the rotation rate of the reference frame relative to an inertial 
framing. Since (3)-(5) are independent of 12, the standard K-e model predicts the some 
results for all rotation rates and, hence, does not distinguish between turbulent shear 
flow in an inertial frame and rotating turbulent shear flow. Speziale and Mac Giolla 
Mhuiris [11] recently showed that the K-e model has the following equilibrium solution 
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for rotating shear flow: 


= (M~ = -j(<W. (i>ia)~ = 0 

(622)00 = 0, (623)00 = 0, (633)00 = 0 

(¥)-(*)' 


(7) 

( 8 ) 

(9) 


where a = (C e 2 — l)/(C e i — 1) and (•)<» denotes the equilibrium value obtained in the 
limit as t — > 00 . These equilibrium values are universal, i.e., are completely independent 
of the initial conditions, the shear rate, and the rotation rate. It was also shown in [11] 
that the long time solutions (corresponding to t* = St 1) for the kinetic energy and 
dissipation rate in the K-e model grow exponentially: 


K ~ exp 


e ~ exp 


l£±( a -l)f 

a 



( 10 ) 


( 11 ) 


Hence, the K-e model predicts the following physical picture for rotating shear flow: the 
turbulent kinetic energy and dissipation rate grow exponentially in time at a compa- 
rable rate; the anisotropy tensor 6# and shear parameter SK/e approach a universal 
equilibrium. While this characterization is qualitatively correct for pure shear flow (see 
Tavoularis and Corrsin [12]), it is quite incorrect for most values of Q/S in rotating shear 
flow. Linear stability analyses and numerical simulations of the Navier-Stokes equations 
indicate that for values of the rotation rate that are discernibly outside of the range 
0 < Cl/S < 0.5, the flow undergoes a restabilization wherein K and t — ► 0 as t — * 00 
(see Bardina, Ferziger, and Reynolds [13] and Bertoglio [14]). 

It will now be demonstrated that, unlike the commonly used two-equation models, 
second-order closures are able to describe the stabilizing or destabilizing effect of rota- 
tions on turbulent shear flow. Speziale and Mac Giolla Mhuiris [11] recently considered 
a fairly general class of second-order closure models of the form: 


■ d^v ' 2 

Tij = ~l" Tjk^mki^m') 

■ r ± r 9% _ r t_ 

el K 13 dx , c2 K 


( 12 ) 


(13) 


for any rotating homogeneous turbulence. In (12)-(13), Ily denotes the pressure-strain 
correlation which is assumed to be of the general form 


Htj = n tJ (r tJ , 


(14) 
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and G e \ and C e2 are either constants or functions of the invariants of bij. This class 
of second-order closures encompasses a wide variety of models including the simplified 
form of the Launder, Reece, and Rodi model for which 


n* = -f —KSij) — c 2 


7~ik I" £mkj*lr 


fdUi 

I -f- £mki^7] 




(15) 


and Ci = 1.8, C 2 = 0.6, C c \ = 1.44 and C e2 = 1.92. It was shown by Speziale and Mac 
Giolla Mhuiris [11] that this class of second-order closure models has two-equilibrium 
solutions for rotating shear flow: one where 



(16) 


which exists for all 0,/S and one where 



7i +72 




(17) 


which exists for a small intermediate band of fi/5 which can range from —0.1 < Cl/S < 
0.6 (here, 70,71, and 72 are directly related to the constants of the model). The former 
equilibrium solution for which (e/SK ) 0 0 = 0, predominantly is connected with solutions 
wherein K and e undergo a power law decay in time; the latter equilibrium solution 
(17), where (e/SK) O0 is nonzero, is connected with solutions where K and e grow expo- 
nentially in time at the same rate. In this intermediate band of tl/S, these two solutions 
exchange stabilities in a fashion that qualitatively mimics the shear instability with its 
exponential time growth of disturbance kinetic energy. 

In Figure 4(a), a bifurcation diagram is shown for the Launder, Reece, and Rodi 
model. This bifurcation structure qualitatively mimics the stabilizing or destabilizing 
effects of rotations on turbulent shear flow as discussed above. In stark contrast to the 
bifurcation that is properly predicted by the second-order closure, the equilibrium dia- 
gram for the K-e model shows the erroneous prediction of a universal value for ( e/SK ) OQ 
which is completely independent of 0,/S (see Figure 4(b)). As mentioned earlier, this 
universal equilibrium solution for the K-e model corresponds to an unstable flow wherein 
K and e grow exponentially in time. 

In addition to yielding a superior qualitative description of the equilibrium structure 
of rotating turbulent shear flows, the quantitative values of the equilibrium states for pure 
shear flow predicted by the second-order closures are also substantially better than those 
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obtained from the commonly used two-equation models. To illustrate this superiority 
of the second-order closures, the equilibrium values for bij and SK/e obtained from the 
Launder, Reece, and Rodi model and the K-e model are compared in Table 1 with the 
experimental results of Tavoularis and Corrsin [12] for homogeneous turbulent shear 
flow. 

In Figures 5(a)-(c), the time evolution of the turbulent kinetic energy predicted by the 
Launder, Reece, and Rodi model and the K-e model are compared with results from the 
large-eddy simulations of Bardina, Ferziger, and Reynolds [13] for three rotation rates: 
Cl/S = 0, Cl/S = 0.25, and Cl/S = —0.5. A direct comparison of Figure 5(a) and Figure 
5(b) with the large-eddy simulations shown in Figure 5(c), graphically demonstrates 
the superior capability of second-order closure models in predicting the stabilizing or 
destabilizing effect of rotations on shear flow. 


3. Needed Modeling Improvements in Second- Order 
Closures 


While the author is in full agreement with the main points of the Launder position 
paper concerning the superior capabilities of second-order closure models, it must be 
cautioned that these models have not yet matured to the point where reliable predictions 
can be made for a variety of complex turbulent flows. Several areas where improvements 
are needed (some of which were pointed out by Professor Launder), will be discussed in 
more detail in this section. 

Second-order closure models are based on the Reynolds stress transport equation 
which takes the exact form 


dTjj 

dt 


, _ dtij dvj dv { dCij k 

+ Vk dx k ~ Tlk dx k Tjk dx k + dx k 


- n ij + e^ + uv\j 


where 


— it iUjU k + pUiSjk d" pUj8i k 


(18) 

(19) 

( 20 ) 

dx k dx k 

are the third-order diffusion correlation, the pressure-strain correlation, and the dissi- 
pation rate correlation, respectively (u is the kinematic viscosity of the fluid). In order 
for closure to be achieved at this “second moment” level (which forms the raison d’etre 
of second-order modeling), models must be developed wherein the higher-order corre- 
lations Cij k ,Kij and £ij are taken to be functionals of the Reynolds stress mean 


_ du 1 . 

_ 2 TS^ 
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velocity gradients &Ui/dxj and some length scale of turbulence A. The Reynolds stress 
is decomposed into isotropic and deviatoric parts as follows 

Tij = -\k6„ - 2Kbij ( 22 ) 

and the turbulence length scale A is usually assumed to be of the form 

A = C*—~ (23) 

e 

where C* is a dimensionless constant and e = is the turbulent dissipation rate. 
Hence, consistent with the use of (23), the higher-order correlations can be taken to be 
functionals of bij,dvi/dxj, K, and e instead. 

Typically, the third-order diffusion correlation is modeled by a gradient transport 
hypothesis wherein it is assumed that is of the general form 

C iik =C, jUm (b,K,e)^. (24) 

Motivated by analyses based on homogeneous turbulence [15], virtually all of the com- 
monly used models for the pressure-strain correlation are assumed to be of the form 

^+^ =£ ^ (b)+ ^“ (b) S (25) 


where the first term on the right-hand-side of (25) represents the slow pressure strain 
while the second term represents the rapid pressure strain. Typically, the turbulence 
dissipation correlation is assumed to be of the general form 


£ij — 2 ® d* fa^ij 


(26) 


where f, is taken to be dimensionless constant or a function of the invariants of bij. For 
high Reynolds number turbulent flows that are sufficiently far from solid boundaries, f„ 
is taken to be zero. In order to achieve closure of the Reynolds stress transport equation, 
(24)-(26) must be supplemented with a modeled transport equation for the turbulent 
dissipation rate which is of the general form 

d £ d £ » . . 

fo+v k —~vV 2 e + V c -$ e + V e (27) 

where V e , $ e and V e represent the production, dissipation, and turbulent diffusion of e. 

Most of the existing second-order closure models can be constructed by expanding 
the unknown tensor coefficients on the right-hand-sides of (24)-(25) in a Taylor series in 
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bij (subject to the symmetry properties of C#* and 11#). The older models are actually 
first-order Taylor expansions in bij] for example in the Launder, Reece, and Rodi model 

M= - „ 

OTjk 


C ijk = -C,-(r u 


^Tik d Tij v 

dxi Tj7 dxi Tkl dxi 


n# = -C 1 eb ij + C 2 KS ij + C 3 K(b ik S jk 

+b jk s ik - ib^SnJij) + C 4 K(b ik W jk + b jk W ik ) 


(28) 

(29) 


where 5# = K&uJdxj + dvj/dxj ) and W# = \ (dvi/dxj - dvj/dxi ) are the mean rate of 
strain tensor and vorticity tensor, respectively. In the simplified version of the Launder, 
Reece, and Rodi model (which is now referred to as the “Basic Model” by Launder and 
his co-workers), C\ = 1.8, C 2 = 0.8, C 3 = 1.2, C 4 = 1.2, and C, = 0.11. 

The modeled terms in the dissipation rate transport equation (27) are typically based 
on the assumption that the production (or dissipation) of the turbulent dissipation 
is proportional to the production (or dissipation) of the turbulent kinetic energy. A 
gradient transport hypothesis is typically invoked for the turbulent diffusion term on the 
right-hand-side of (27). These assumptions give rise to a modeled transport equation 
for the turbulent dissipation rate that takes the general form 

dvi 


de _ de 

dt Vt dxi 


- C * K r 'i dx . 


-cj= - 


K 


(30) 


d_ 

dxi 


. _ K de . 

(c, 7 r «&; ) + " v '■ 


In the Launder, Reece, and Rodi model, C e 1 , C e2 , and C e are taken to be constants which 
assume the values of 1.44, 1.92, and 0.15, respectively. Some more recent models have 
taken C e i and C e2 to be functions of some subset of the invariants of 6# and &Uj/dxj] 
these newer models will be discussed in more depth in the next section. 

Now, with the aid of this background material, the four active areas of research for 
the development of improved second-order closure models that were mentioned in the 
introductory remarks can be elaborated on. These areas are as follows: 

(i) The development of improved models for the pressure-strain correlation of turbu- 
lence which account for nonlinear anisotropic effects. Since for most flows of engineering 
interest, ||b|| ~ 0.2, the use of a first-order Taylor expansion in b for A# and is 

highly questionable. In fact, nonlinear terms in the model for A# axe needed to predict 
the curved trajectories that occur in the phase space of the return to isotropy problem 
as shown in Figure 6 (see Choi and Lumley [9] and Sarkar and Speziale [17]). At least 
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a quadratic nonlinearity in the model for Mijki is needed in order to satisfy the con- 
straint of Material Frame Indifference (MFI) in the limit of two-dimensional turbulence 
(Speziale [18] and Haworth and Pope [19]). This MFI constraint is the mathematical 
embodiment of the well-known result that two-dimensional disturbances evolve iden- 
tically in both a rotating frame and an inertial frame.* In addition, recent work on 
rotating turbulent shear flows (Speziale, Sarkar, and Gatski [20]) and Rapid Distortion 
Theory (Reynolds [21]) have suggested the possible need for new terms in (25) that are 
nonlinear in the mean velocity gradients &Di/dxj. 

(ii) The development of non-gradient transport models for the turbulent diffusion 
terms such as Cijk need to be considered seriously. It is well-known that turbulent flows 
do not have a clear cut separation of scales; the largest eddies (which contain a significant 
portion of the turbulent kinetic energy) are of a comparable size to the geometrical scale 
of the flow. Consequently, one would expect a gradient transport hypothesis (which is 
rigorously derived as a first-order expansion in the ratio of fluctuating to mean length 
scales) to only constitute a crude approximation. Many difficulties in the prediction 
of turbulent mixing layers could be tied to this deficiency in the modeling of turbulent 
diffusion by means of gradient transport. 

(iii) The development of asymptotically consistent near wall corrections to the tur- 
bulence models for 11; j, £ij and C^k that are geometry-independent and do not have any 
ad hoc damping functions. Most of the commonly used corrections are either asymp- 
totically inconsistent, geometry-dependent through an artificially imposed dependence 
on the unit normal to the wall, or otherwise ad hoc through the use of wall damping 
functions based on the turbulence Reynolds number or the distance from the wall (see 
Hanjalic and Launder [22]). Such empiricisms do not allow for the reliable prediction 
of wall transport properties (e.g., skin friction and heat transfer coefficients) that are 
extremely important in aerodynamic applications. 

(iv) The development of improved modeled transport equations for the turbulence 
length scale is an issue of utmost importance. Most of the commonly used models have 
a scalar turbulence length scale A based on the dissipation (A <x iff /e). There is the 
obvious objection that the construction of a turbulence macro-scale based on small-scale 
(one-point) information is conceptually wrong. Furthermore, this rather simplified defi- 
nition of length scale contains no directional information. Although attempts have been 
made to develop a length scale equation based on an integral of the two-point velocity 
correlation tensor (Wolfshtein [23] and Donaldson and Sandri [24]), it can be shown that 

* Consequently, the models that violate this constraint cannot be used in the analysis of geostrophic 
turbulence. 
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these specific models are equivalent to the standardly used dissipation rate transport 
model in the limit of homogeneous turbulence - a simplified case for which the modeled 
dissipation rate transport equation is already deficient. Fundamentally new research is 
needed on the development of modeled transport equations for some appropriate choice 
of integral length scales that contain the required directional information. 

Finally, in regard to these four points, it should be mentioned that some significant 
progress has been made in the development of improved models for the rapid pressure- 
strain correlation by means of realizability (Lumley [25], and Shih and Lumley [26]), and 
the invariance considerations discussed in point (i). However, in the opinion of the au- 
thor, very little progress, if any, has been made since the early 1970’s in the development 
of better models for the turbulence length scale and diffusion effects. An examination 
of the modeled dissipation rate transport equation as a basis for the turbulence length 
scale will be discussed in the next section. 


4. The Modeled Dissipation Rate Transport Equa- 
tion 


Now, the strengths and weaknesses of the commonly used modeled dissipation rate 
transport equation will be discussed. Furthermore, an attempt will be made to demon- 
strate at what level of approximation this model is derivable from the two-point corre- 
lation tensor which more properly contains information about the turbulent macroscale. 
Finally, an argument will be put forth as to why attempts at the development of im- 
proved modeled dissipation rate transport equations have failed during the past decade. 

In order not to cloud the issue with the added difficulties that are associated with the 
integration of turbulence models to a solid boundary, the more simplified case of homo- 
geneous turbulence will be considered. For any homogeneous turbulence, the standardly 
used version of the modeled dissipation rate transport equation (30) reduces to 


. e dvi 

C ‘' K Tii dx i 



(31) 


The transport equation for the turbulent kinetic energy has the exact form 


ir_ dv< 
K ~ Tii dxj 


(32) 


and is obtained from a contraction of (18). Equations (31)-(32) can be combined to 
yield a transport equation for the reciprocal turbulent time scale e / K which is given by 


d_ 

dt 


G?) = 2(1 " + (1 " c ’ 2) G?) 


(33) 
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Provided that C e \ > 1, this equation has two equilibrium solutions in the limit as t — * oo: 

(*L = ° < 34 > 




dvi 

'Ifoj 


(35) 


which are obtained by setting the time derivative on the left-hand-side of (33) to zero. 
C c \ = 1 constitutes a bifurcation point of equation (33); for C e i < 1, the only realizable 
fixed point of (33) is (e/if) TO = 0. The zero fixed point (34) predominantly corresponds 
to solutions for if and e that undergo a power law decay in time [11] (i.e., e and if — > 0 
as t — * oo). This branch of solutions allows for the prediction of isotropic decay as well 
as the flow restabilization that occurs at certain rotation rates in rotating plane shear 
and plane strain turbulence. Furthermore, since from (33) we have 

d 


dt \K) 


0 


(36) 


when e/if = 0, and from (31) we have 


^ = 0 
dt 


(37) 


when e = 0, we conclude that if e(0) > 0 and if(0) > 0 then e(t) > 0 and K{t) > 0 
for all later times t. Hence, the standardly used modeled e-transport equation guarantees 
realizability with respect to if and e by virtue of the fact that e/if = 0 is, in dynamical 
systems terms, an “invariant plane.” 

The non-zero fixed point (35) is associated with solutions for K and e which grow 
exponentially in time [11]. More precisely, for the non-zero fixed point (35), it can be 
shown that 

if ~ exp (At), e ~ exp(At) (38) 

for \t 1, where 

*-ww-£+Gf)j- w 

This allows for the prediction of a structural equilibrium in homogeneous turbulent shear 
flow wherein SK/e achieves an equilibrium value that is independent of both the initial 
conditions and the shear rate. Such an equilibrium for turbulent shear flows has been 
observed experimentally by Tavoularis and Corrsin [12] for weak to moderately strong 
shear rates. 

It is thus clear that the two fixed points (34)-(35) of the commonly used e-transport 
equation have certain properties that are crucial to the proper description of homoge- 
neous turbulent flows. Now, it will be shown that recently proposed alterations to this 
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modeled e-transport equation destroyed one or the other of these key fixed points and, 
hence, were doomed to failure. 

For example, based on the desire to account for rotational strains, Pope [27] proposed 
an alteration to the e- transport equation wherein a term of the form 

c, 3 ^-'s„w jk w ki 

(where (7*3 is a constant) was added to the right-hand-side of (30). This eliminates the 
(e/if )oo = 0 fixed point, for nonzero *S’ tJ W jk W & , which can cause problems with realiz- 
ability and can eliminate the ability to predict flow restabilizations in three-dimensional 
rotating shear flows. Hanjalic and Launder [28] proposed a modification wherein the 
term 

-CaKWijWij 

was added to the right-hand-side of (30) (where C e 3 is a different constant). It is a simple 
matter to show that, for any nonzero Wij, this alteration also eliminates the fixed point 
(e/ K) oo = 0. Consequently, this model can have problems with realizability and predicts 
that the turbulent kinetic energy and dissipation rate oscillates (for all times t > 0) in 
rotating isotropic turbulence - a prediction that is in substantial contradiction to the 
results of physical and numerical experiments which suggest a monatonic power law 
decay. Most recently, Launder and his co-workers proposed a modified version of the 
e-transport equation for which 

C. l = l,C a = C a {lI,III). (40) 

In this case, the choice of C e \ — 1 eliminates the nonzero fixed point (35) (the only fixed 
point is (e/if)oo = 0). Consequently, this model will not predict an exponential time 
growth of K and e in homogeneous shear flow (and SK/e will not approach a universal 
equilibrium) in apparent contradiction to physical and numerical experiments. Similar 
problems occur in rotating turbulent shear flows with a modification proposed recently 
by Bardina, Ferziger, and Rogallo [29]. 

It has been demonstrated that the standard form of the modeled dissipation rate 
transport equation has several crucial properties that have been destroyed by virtually 
every attempt to modify it. To understand this point more deeply, it would be helpful 
to see what approximations are necessary to derive the modeled e-transport equation 
(31) from an analysis of the dynamics of the two-point correlation tensor 

Rij = Ui(x)uj(x + r). (41) 
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A study conducted not too long ago by Donaldson and Sandri [24] can shed some light 
on this issue. They introduced a tensor length scale A ij defined by the solid angle 
integration 

A modeled transport equation for A ij was obtained as follows: (a) the exact transport 
equation for Rj was closed by assuming that the higher-order unknown correlations 
could be expanded in a Taylor series in (subject to the appropriate symmetry prop- 
erties) where only first-order terms were kept, and (b) the resulting modeled transport 
equation for Rj was multiplied by l/47rr 2 and integrated over all of r space. The de- 
rived transport equation for A,j obtained by this method is as follows for homogeneous 
turbulence [24]: 

U K A*> - -* a *§-* a *£ 

+-/2v a ^-K ij - - A Si,) (43) 

-2V2 bK 5 Sij 

where A = |A,-,; v C3 and b are constants. The modeled Reynolds stress transport 
equation that is solved in conjunction with (43) is obtained from a simple contraction 
of the modeled transport equation for R{j\ this equation takes the form [24] 

drij dvj dvi r-K , 2 r _ . . 

dt - ^ rit gx k ~ Tik dx k ~'^J {T ‘ i + 3 KS ‘ i) 

'(«) 

+ 3 ^ b ~A Sij - 

From a direct comparison of the contraction of (44) with (32), it follows that 

e = 2y/2 b—r (45) 

A 

for this tensor length scale model in homogeneous turbulence. By taking the time 
derivative of (45), the transport equation 

„ k ,,3 Ki • K* .. 

1 e &Si y/2 e 2 dvi , . 

= 2K t< + 6 (46) 


_/3 

l 2 2 b’K 
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is obtained. If an anisotropy tensor is defined for the tensor length scale in the following 
manner 

A ... _ 1 A... - 

(47) 


t(a) _ 3 

* “ Akk 

it can be shown that b\^ is a solution of a transport equation that is of exactly the same 
form as that for bij. Hence, if we consider a homogeneous turbulence that evolves from 
an initially isotropic state, then 

= bij (48) 


and 


A 


K 5_ dvj 

dxj 


'XJ 


(49) 


Then the dissipation rate transport equation (46) corresponding to this tensor length 
scale model simplifies to 


. e ovi o, u C2 .w 
e = 1 - 5 ^ r *i— + 


dVj 




Vo 

3 b'K 


which is of the same form as the standardly used equation (31) with C e i 
Ce2 = (§ + v ca /2b) » 1.80. 


(50) 
1.5 and 


It has thus been shown that the standard form of the modeled e-transport equation 
can be derived from the two-point correlation tensor by making approximations that are 
comparable to those made in the derivation of the modeled Reynolds stress transport 
equation by Launder, Reece, and Rodi [16]. Consequently, it seems to be somewhat 
questionable to argue (as many have done) that the standard form of the e-transport 
equation is the weak link in the commonly used second-order closures. There is no doubt 
that improvements in the specification of the turbulence length scale are needed, but 
these are more likely to come from general analyses of some set of modeled transport 
equations for the integral length scales which incorporate directional information. The 
kind of ad hoc adjustments in the modeled dissipation rate transport equation that have 
been considered during the past decade appear to be counterproductive. 


5. Concluding Remarks 

Projected advances in computer capacity make it highly unlikely, in the opinion of 
the author, that the complex turbulent flows of engineering interest will be solved on a 
routine basis by direct or large-eddy simulations for at least several decades to come. 
Furthermore, substantial theoretical difficulties with two-point closures for complex tur- 
bulent flows that are strongly inhomogeneous make their application to problems of 
engineering interest equally unlikely for the foreseeable future. It appears that Reynolds 
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stress models are likely to remain the method of choice for the solution of the turbulence 
problems of technological interest for at least the next few decades. Direct numerical 
simulations of the Navier-Stokes equations will continue to be restricted to less geomet- 
rically complex turbulent flows, at lower Reynolds numbers, where they will be used to 
gain a better insight into the basic physics of turbulence. 

Among the variety of existing Reynolds stress closures, it is the opinion of the author 
that second-order closure models constitute, by far, the most promising approach. The 
predictive capabilities of second-order closures have been enhanced somewhat over the 
past decade due to significant modeling improvements in the pressure-strain correlation. 
However, new prescriptions for the turbulent length scale (which incorporate some di- 
rectional and two-point information) as well as new methods for integration to a solid 
boundary are direly needed before more reliable models can be obtained. In fact, these 
issues are of such overriding importance for wall-bounded turbulent flows that Reynolds 
stress model predictions can be degraded to the point where they are no better than 
those of the K-e model - a state of affairs that has contributed to some of the misleading 
critical evaluations of second-order closures that have been published in the recent past. 

It must be remembered that second-order closure models are one-point closures and 
thus can never yield accurate quantitative predictions for a wide variety of turbulent flows 
where the energy spectrum can change drastically. Nonetheless, with the implementation 
of the improvements discussed in this paper, it is quite possible that a new generation of 
second-order closures can be developed that will provide acceptable engineering answers 
for a significant range of turbulent flows that are of technological interest. 
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Equilibrium 

Values 

Standard 
K — e Model 

Launder, Reece 
&: Rodi Model 

Experiments 

(M- 

0 

0.193 

0.201 

(^ 22)00 

0 

-0.096 

-0.147 

(^ 12)00 

-0.217 

-0.185 

-0.150 

(SK/e)oo 

4.82 

5.65 

6.08 


Table 1. Comparison of the predictions of the standard K — e model and the Launder, 
Reece, and Rodi model with the experiments of Tavoularis and Corrsin [12] 
on homogeneous shear flow. 


17 




O Experimental data of Choi and Lumley [9] 
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J 2. Temporal decay of the second invariant of the anisotropy tensor: Comparison 
of the Rotta model with the relaxation from plane strain experiment of Choi 
and Lumley [9]. 







(a) K—t Model 


(t>) LRR Model 


0/S = 0, 0.25 6c -0.50 


n/s =-o.5o 


0/s = 0.25 


0/S = 0 


0/s =-0.50 


Time evolution of the turbulent kinetic energy for homogeneous turbulent 
shear flow in a rotating frame: (a) K-e model, (b) Launder, Reece and Rodi 
model, (c) Large-eddy simulations of Bardina, Ferziger and Reynolds [13]. 








Figure 6. Phase space diagram of the relaxation from plane strain experiment of Choi 
and Lumley [ 9 ] (£ = III 1 / 3 , 77 = II 1 / 2 ). 
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